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Abstract 

In this paper, a new type of multi-level correction scheme is proposed for 
solving eigenvalue problems by finite element method. With this new scheme, 
the accuracy of eigenpair approximations can be improved after each correc- 
rS^ • tion step which only needs to solve a source problem on finer finite element 

d . space and an eigenvalue problem on the coarsest finite element space. This 

correction scheme can improve the efficiency of solving eigenvalue problems 
by finite element method. 
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C^ I ment method, multi-space, multi-grid. 
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2 1 Introduction 

The purpose of this paper is to propose a new type of multi-level correction scheme 
^ . based on finite element discretization to solve eigenvalue problems. The two-grid 

method for solving eigenvalue problems has been proposed and analyzed by Xu and 
Zhou in [21J. The idea of the two-grid comes from P^ I2U] for nonsymmetric or 
indefinite problems and nonlinear elliptic equations. Since then, there have existed 
many numerical methods for solving eigenvalue problems based on the idea of two- 
grid method ([DEI [IT]). 
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In this paper, we present a new type of multi-level correction scheme for solving 
eigenvalue problems. With the new proposed method solving eigenvalue problem will 
not be much more difficult than the solution of the corresponding source problem. 
Our method is some type of operator iterative method (llHETlEBj). The correction 
method for eigenvalue problem in this paper is based on a series of finite element 
spaces with different approximation properties which are related to the multilevel 
method ([IB]). 

The standard Galerkin finite element method for eigenvalue problem has been 
extensively investigated, e.g. Babuska and Osborn [21 [3], Chatelin [5] and references 
cited therein. Here we adopt some basic results in these papers in our analysis. The 
finite element method for eigenvalue problem has been developed well and many 
high efficient methods have also been proposed and analyzed for different types of 
eigenvalue problems ([HlHIIOllISlIIllIiailllEI]). 

The corresponding error estimates of the type of multi-level correction scheme 
which is introduced here will be analyzed. Based on the analysis, the new method 
can improve the convergence rate of the eigenpair approximations after each cor- 
rection step. The multi-level correction procedure can be described as follows: (1) 
solve the eigenvalue problem in the coarsest finite element space; (2) solve an addi- 
tional source problem in an augmented space using the previous obtained eigenvalue 
multiplying the corresponding eigenfunction as the load vector; (3) solve eigenvalue 
problem again on the finite element space which is constructed by combining the 
coarsest finite element space with the obtained eigenfunction approximation in step 
(2). Then go to step (2) for next loop. 

Similarly to [21], in order to describe our method clearly, we give the following 
simple Laplace eigenvalue problem to illustrate the main idea in this paper with 
multi-grid implementation way (see section 5). 

Find(A, u) such that 




Xu, in Q, 
0, 

1, 



0, on on, (1.1^ 



where fi C 7^^ is a bounded domain with Lipschitz boundary dQ and A denote the 
Laplace operator. 

Let Vh denote the coarsest linear finite element defined on the coarsest mesh Th- 
Additionally, we also need to construct a series of finite element spaces Vh^, V^g, 
■ ■ ■ ) '^h„ which are defined on the corresponding series of meshes 7/^. (fc = 2, 3, ■ ■ ■ n) 
such that Vh C Vh^ C ■ ■ ■ C Vh„ ([HE]). Our multi- level correction algorithm to 
obtain the approximation of the eigenpair can be defined as follows (see section 3 
and section 4): 

1. Solve an eigenvalue problem in the coarsest space Vh'- 



Find {\h,uh) E TZ x Vh such that llwi^-llo = 1 and 



Jii Jq 



2. Set hi = H and Do k = 1,- ■ ■ , n — 2 

• Solve the following auxiliary source problem: 
Find Mhfc+i e Vhk+i such that 



Ja Jq 

• Define a new finite element space VH,hk+i = Vh + spanJM/ij.^^ } and solve 
the following eigenvalue problem: 

Find {\h^^^,Uh^^J eTZx VH,hk+i such that ||Mfe,+J|o = 1 and 

/ ^Uh^+yVH,h,+,dn = Xh^^^ / Uh^^^VH,h^+,dn, yVHM+i e Vf/,hfe+i- 

Jn Jn 



end Do 

3. Solve the following auxiliary source problem: 
Find Uh„ € Vh such that 



Vuh^Vh^dn = Xh„_, / Uh„_,Vh„dn, Wvh„ e Vh„. 
n Jn 

Then compute the Rayleigh quotient 

IIV^!, Iln 

M il'n W U 



A 



'" "uhM 



If, for example, A^ is the first eigenvalue of the problem at the first step and 1] is a 
convex domain, then we can establish the following results (see section 3 and section 
4 for details) 



n n 



r2(n-fc) 
fc=l ■ ' k=l 



These two estimates means that we can obtain asymptotic optimal errors by taking 
H = 01^ andhk = H'' {k = 2,--- ,n-l). 

In this method, we replace solving eigenvalue problem on the finest finite element 
space by solving a series of boundary value problems in the corresponding series 
of finite element spaces and a series of eigenvalue problems in the coarsest finite 



element space. As we know, there exists multigrid method for solving boundary 
value problems efficiently. So this correction method can improve the efficiency of 
solving eigenvalue problems. 

An outline of the paper goes as follows. In Section 2, we introduce finite element 
method for eigenvalue problem and the corresponding error estimates. A type of one 
correction step is given in section 3. In Section 4, we propose a type of multi-level 
correction algorithm for solving eigenvalue problem by finite element methods. In 
section 5, two numerical examples are presented to validate our theoretical analysis 
and some concluding remarks are given in the last section. 

2 Discretization by finite element method 

In this section, we introduce some notation and error estimates of the finite element 
approximation for eigenvalue problems. In this paper, the letter C (with or without 
subscripts) denotes a generic positive constant which may be different at different 
occurrences. For convenience, the symbols <, > and ~ will be used in this paper. 
That xi < yi,X2 > y2 and X3 ~ 1/3, mean that xi < Ciyi, X2 > C2?/2 and C3X3 < 
2/3 < C'sa^s for some constants Ci, C2, C3 and C3 that are independent of mesh sizes. 

Let (y, II ■ II) be a real Hilbert space with inner product (-,■) and norm || ■ ||, 
respectively. Let a(-, ■), &(■,■) be two symmetric bilinear forms on X ^ X satisfying 

aiw^v) < ||i(;|| ||f II, V-u; G l^ and Vi) G V, (2.1) 

ll^lP ~ ci{w,w), Vti? G y and < &(w, tf), Vto G \^ and if 7^ 0. (2.2) 

From 02. ip and (12. 2p . we know that || ■ ||a := ci{-^ 'Y^'^ and || ■ || are two equivalent 
norms on V . We assume that the norm || ■ || is relatively compact with respect to 
the norm || • \\b := b{-, ■)^/^. We shall use a(-, ■) and || ■ ||a as the inner product and 
norm on V in the rest of this paper. 

Set 

W := the completion of V with respect to || ■ \\b. 

Thus VF is a Hilbert space with the inner product b{-, •) and compactly imbedded in 
V. Construct a "negative space" by V = the dual of V with a norm || ■ ||_a given 
by 

||i(;||_a = sup b{w,v). (2.3) 

VGV,\\v\\a = l 

Then W C V compactly, and for v E V, b{w,v) has a continuous extension to 
w E V such that b{w,v) is continuous on V by Hahn-Banach theorem ([9]). We 
assume that Vh G V is a. family of finite-dimensional spaces that satisfy the following 
assumption: 



For any w ^ V 



lim inf Iliy — v|L = 0. (2.4) 

h-^Ov£Vh 



Let Ph be the finite element projection operator of V onto Vh defined by 

a{w - PhW, v) = 0, Vu; G y and ^v G 14- (2.5) 

Obviously 

\\PhW\\a<\\w\\a, WweV. (2.6) 

For any w & Vjhj (12. 4p we have 

\\w-Phw\\a = o(l), as/i^O. (2.7) 

Define ria{h) as 

Va{h)= sup inf ||T/-i;|U, (2.8) 

/eV,||/IU=i^el4 

where the operator T : V ^ V is defined as 

a{Tf,v) = b{f,v), V/ G V and Vt; G V. (2.9) 

In order to derive the error estimate of eigenpair approximation in negative norm 
II • II _a, we need the following negative norm error estimate of the finite element 
projection operator Ph. 

Lemma 2.1. (J^ Lemma 3.3 and Lemma 3.4]) 

r]a{h) = o(l), as/i^O, (2.10) 

and 

\\w-PhW\\^a < Vaih)\\w - PhW\\a, ^W ^ V. (2.11) 

In our methodology description, we are concerned with the following general eigen- 
value problem: 

Find {X, u) E TZ X V such that b{u, u) = 1 and 

a{u,v) = \b{u,v), yveV. (2.12) 

For the eigenvalue A, there exists the following Rayleigh quotient expression ([2], 

El ED) 

a(u, u) 
A = -^. 2.13 

b[u,u) 
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From [21[S]) we know the eigenvalue problem fl2.12p has an eigenvalue sequence {Xj} : 
< Ai < A2 < ■ ■ ■ < Afc < • ■ ■ , lim Afc = 00, 

fc— )-oo 

and the associated eigenfunctions 

ui,U2,--- ,Uk,--- , 

where b{ui,Uj) = 5ij. In the sequence {A^}, the \j are repeated according to their 
geometric multiplicity. 

Now, let us define the finite element approximations of the problem (12.120 . First 
we generate a shape-regular decomposition of the computing domain Vt C Tl!^ {d = 
2, 3) into triangles or rectangles ioT d = 2 (tetrahedrons or hexahedrons for d = 3). 
The diameter of a cell K & Th is denoted by hx- The mesh diameter h describes 
the maximum diameter of all cells K E Th- Based on the mesh Th, we can construct 
a finite element space denoted hj Vh G V. In order to do multi-level correction 
method, we start the process on the original mesh Th with the mesh size H and the 
original coarsest finite element space Vh defined on the mesh Th- 

Then we can define the approximation of eigenpair {X,u) of fl2.12p by the finite 
element method as: 

Find (Aft, Uh) eTZ x Vh such that b{uh, Uh) = 1 and 

a(uh,Vh) = Xhb{uh,Vh), WvheVh. (2.14) 

From (I2.14p . we can know the following Rayleigh quotient expression for A^ holds 

(lasiizii) 

A. = 'p^. (2.15) 

b[Uh,Uh) 

Similarly, we know from [3l |5] the eigenvalue problem (I2.12p has eigenvalues 
and the corresponding eigenfunctions 

Ul,h, U2,h, ■ ■ ■ , Uk,h, ■ ■ ■ , UNf^h, 

where b{ui^h,Uj^h) = ^ijA < ^,J < Nh {Nh is the dimension of the finite element 
space Vh). 

From the minimum- maximum principle ([21 13]), the following upper bound result 
holds 

Xi < Xi^h, i = l,2,--- ,Nh. 



Let M(Aj) denote the eigenspace corresponding to the eigenvalue Aj which is defined 
by 

M{Xi) = [w G V : w is an eigenvalue of (12.121) corresponding to Aj 

and ||w||6 = 1}. (2.16) 

Then we define 

Sh{>^i) = sup inf \\w - v\\a. (2-17) 

For the eigenpair approximations by finite element method, there exist the fol- 
lowing error estimates. 

Proposition 2.1. (;2, Lemma 3.7, (3.29b)],l3, P- 699] and f^) 

(i) For any eigenfunction approximation Ui^h of i2.14\) {i = 1,2, ■■ ■ ,Nh), there is 
an eigenfunction Ui of /i2.12\) corresponding to Aj such that \\ui\\f, = 1 and 

\\u.-Ui4a < CM^i)- (2.18) 

Furthermore, 

\\Ui - Ui^h\\-a < Cir]a{h)\\Ui - Ui^hWa- (2.19) 

(a) For each eigenvalue, we have 

Xi < Xi,h <Xi + Ci5l{\) (2.20) 

Here and hereafter Ci is some constant depending on i hut independent of the mesh 
size h. 



3 One correction step 

In this section, we present a type of correction step to improve the accuracy of 
the current eigenvalue and eigenfunction approximations. This correction method 
contains solving an auxiliary source problem in the finer finite element space and 
an eigenvalue problem on the coarsest finite element space. For simplicity of no- 
tation, we set (A,m) = (Aj,Mj) {i = 1,2,- •■ , fc, ■ ■ ■ ) and {Xh,Uh) = {Xi,h,Ui^h) {i = 
1, 2, ■ ■ ■ ,Nh) to denote an eigenpair of problem (12.121) and (I2.14p . respectively. 

To analyze our method, we introduce the error expansion of eigenvalue by the 
Rayleigh quotient formula which comes from [2], [3l [161 EI] • 

Theorem 3.1. Assume (A,u) is the true solution of the eigenvalue problem ^2.12\) . 
j^ ip E V. Let us define 



Then we have 

Proof. First from fl2.13p . (13. ip and direct computation, we have 



A-A 



a('0 — «,'?/' — «) + 2a('0, u) — a(u, u) — A&(7/;, ■?/') 
a('0 — u,ip — u) + 2Xb{i/j, u) — Xb{u, u) — Xb{i/j, ip) 
a('0 — u,ip — u) — Xb{i/j — u,ip — u) 



(3.3) 



6(^,^) 

Then we obtain the desired result (13. 2p . D 

Assume we have obtained an eigenpair approximation (A/j^, -u/jj ElZ x Vh^. Now 
we introduce a type of correction step to improve the accuracy of the current eigen- 
pair approximation (A/u,w/iJ. Let V/ij C V^ be a finer finite element space such 
that Vh^ C Vh2- Based on this finer finite element space, we define the following 
correction step. 

Algorithm 3.1. One Correction Step 

1. Define the following auxiliary source problem: 
Find Uh2 e Vh2 such that 

aiuh2,Vh2) = XhMuh^,Vh2), ^Vh^eVh^- (3.4) 

Solve this equation to obtain a new eigenfunction approximation Uh2 G V/ij . 

2. Define a new finite element space Vu^hi = ^h + spanjw/jj} ^'^'^ solve the fol- 
lowing eigenvalue problem: 

Find (A/ijjtt/iJ G 7?. x VH,h2 such that b{uh2,Uh2) = 1 and 

a{uh2,VH,h2) = ^h2b{uh2,VH,h2)^ yvH,h2 ^^H,h2- (3-5) 

Define {Xh2,Uh2) = Correction{VH, Xh^,Uh^,Vh2)- 

Theorem 3.2. Assume the current eigenpair approximation {Xh^,Uhj) G 7^ x V/^ 
has the following error estimates 

\\u-Uh^\\a < £:/m(A), (3.6) 



||M-M/iil|-a < r]aiH)\\u - Uh^Wa, (3.7) 

|A-A,J < e^^(A). (3.8) 

r/ien a/^er one correction step, the resultant approximation {Xh2, Uh2) G 7^ x Vh^ has 
the following error estimates 

Wu-Uh^Wa < £h2(A), (3.9) 

Wu-Uh^W-a < r]a{H)\\u - Uh^Wa, (3.10) 

|A-A,,| < 4(A), (3.11) 

where ^/.^(A) := r]a{H)eh,{X) + el^{X) + dh^iX). 

Proof. From problems ^^, (I^TT^ and ([S3D, and ([SSD, (E2D and (jSS]), the following 
estimate holds 

\\uh2 - Ph2u\\l ^ a-{uh2 - Ph2U, Uh2 - Ph2u) = b{Xh^Uh^ - Xu, Uh2 - -P/ia") 

< ||A/iiM/ji - Xu\\_a\\Uh2 - Ph2u\\a 

^ (|A/ii - A|||U;,J|_„ + A||n;,^ - u\\^a)\\Uh2 " ^/i2^IU 

^ i^lM) + Va{H)ehA>))\\Uh2 - Ph2u\\a- 

Then we have 

l|w/.2-n.w|U, < 4(A) + 7;,(//)£,,(A). (3.12) 

Combining (I3.12p and the error estimate of finite element projection 

\\U-Ph^u\\a < Sh2{X), 

we have 

\\uh2-u\\a < 6l^{X)+r]a{H)eh,W + 6h2iX). (3.13) 

Now we come to estimate the eigenpair solution {Xh^.u^.^) of problem fl3.5|) . Based 
on the error estimate theory of eigenvalue problem by finite element method ([2113]), 
the following estimates hold 

\\u-Uh2\\a < sup inf \\w - v\\a < \\u - Uh^Wa, (3.14) 

and 

\\u-Uh,,\\-a < rja{H)\\u - Uh2\\a, (3.15) 

where 

?ja{H) = sup inf \\Tf-v\l<r^,{H). (3.16) 

/GV,||/|U=i^eyH,ft2 

From dSlSD, dSHD, (IXT^ and (IXTBD . we can obtain ([S3]) and dXTU]) . The estimate 
fl3.1ip can be derived by Theorem 13. II and fl3.9|) . D 
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4 Multi-level correction scheme 

In this section, we introduce a type of multi-level correction scheme based on the 
One Correction Step defined in Algorithm 13.11 This type of correction method can 
improve the convergence order after each correction step which is different from the 
two-grid method in [21]. 



Algorithm 4.1. Multi-level Correction Scheme 

1. Construct a coarse finite element space Vh and solve the following eigenvalue 
problem: 

Find {Xh, uh) &Ti xVh such that b{uH, uh) = 1 and 

a{uH,VH) = \Hb{uH,VH), VVi/ G Vff. (4.1) 

2. Set hi = H and construct a series of finer finite element spaces Vh^-, ■ ■ ■ , Vh„ 
such that r]aiH) > 4, (A) > 5h,iX) > ■ ■ ■ > (5h„(A). 

3. Dok = 0,l,--- ,n-2 

Obtain a new eigenpair approximation (Xhk+n'^hk+i) G 7^ x V^.^.^^ by a correc- 
tion step 

(V+n«fe.+i) = Correction{VH, XH,Uh^,Vh,^,). (4.2) 

end Do 

4- Solve the following source problem: 
Find Uh^ G Vh„ such that 

a{uh„,VhJ = \h^-Avh^-i^Vh^)^ Vt;,,„Gyft„. (4.3) 

Then compute the Rayleigh quotient of Uh„ 

^^r. - 77 7- 14.4J 

Finally, we obtain an eigenpair approximation {Xh„,Uh^) eTZx Vh„. 

Theorem 4.1. After implementing Algorithm \4.1\ the resultant eigenpair approxi- 
mation {Xh„,Uh„) has the following error estimate 

\\uh„-u\\a < £/i„(A), (4.5) 

|A.„-A| < 4 J A), (4.6) 

n 

where e^M = Z ^a(i/)"-'4,(A). 
fc=i 
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Proof. From T]a{H) > (5/ij(A) > (5/i2(A) > ■ ■ ■ > Sh„{\) and Theorem 13.21 we have 

^kU>^) < VaiH)E,,iX) + S,,^,iX), iorl<k<n-2. (4.7) 

Then by recursive relation, we can obtain 

e.„_,(A) < 77,(//)£,„_,(A) + V.-.(A) 

n— 1 

< 5^r/„(/7)"-'=-i<5,,(A). (4.8) 

fc=i 

Based on the proof in Theorem 13.21 and (14.81) , the final eigenfunction approximation 
Uh„ has the error estimate 

\\uh„-u\\a < £l„_,{X) + Va{H)eh„^,{X) + 5h„{X) 

n 

< Y.'^a{HT''5u,{X). (4.9) 

fc=i 

This is the estimate (14.51) . From Theorem 13. II and (14.91) . we can obtain the estimate 
(US}. D 

5 The application to second order elliptic eigen- 
value problem 

In this section, for example, the multi-level correction method presented in this 
paper is applied to the second order elliptic eigenvalue problem. We also discuss 
two possible ways to implement the multi-level correction Algorithm 14.11 The first 
way is the "two-grid method" of Xu and Zhou introduced and studied in [2T] . The 
second one proposed and studied by Andreev and Racheva in [H [T7] uses the same 
mesh but higher order finite elements. 

In (I2.12p . the second order elliptic eigenvalue problem can be defined by 



a(w, v) = Vm ■ AVvdfl, b{u, v) = puvdVt, 
Jn Jq 

where Q C TZ'^ {d = 2,3) is a bounded domain, A G (W^'°°{Q)) a uniformly 
positive definite matrix on Q and p G W^'°°{Q) is a uniformly positive function 
on Q. We pose Dirichlet boundary condition to the problem and it means here 
V = Hq{Q) and W = L'^{Q). In order to use the finite element discretization 
method, we employ the meshes defined in section 3. 

Here, we introduce two ways to implement the multi-level correction Algorithm 
14.11 The first way uses finer meshes to construct the series of finite element spaces. 

11 



The advantage of this approach is that it uses the same finite element and does not 
require higher regularity of the exact eigenfunctions ([IZ|)- The second way is based 
on the same finite element mesh but using higher order finite elements. In order to 
improve the convergence order, the higher regularity of the exact eigenfunctions is 
required. 

Let us discuss the methods to construct the series of finite element spaces Vh^{k = 
2, 3, ■ ■ ■ , n) for implementing the multi- level correction method. 

Way 1. ("Multi-grid method"): In this case, Vh^ (/c = 2, 3, ■ ■ ■ , n) is the same type 
of finite element as Vh on the finer mesh 7/^. with smaller mesh size hk . Here Th^ is 
a finer mesh of Vt which can be generated by the refinement just as in the multigrid 
method ([21]) from 7/n._i such that hk = ria{H)hk-i. Assume the computing domain 
n is a convex domain. Then rjaiH) = 0{H) and 5h^ = 0{hk) = 0{H^) {k = 
1, 2, ■ ■ ■ , n), and we can obtain the following error estimate for {Xh^,Uh„) 

|A-A,J < VaiHy--'5liX) = OiH'-), (5.1) 

h-unja < Va{Hr-'6H{X) = 0{H-). (5.2) 

From the error estimates above, we can find that the multi-level correction scheme 
can obtain the accuracy as same as solving the eigenvalue problem on the finest 
mesh Th„- This improvement costs solving the source problems on the finer finite 
element spaces Vh^. (A; = 2,3, ■ ■ ■ , n) and the eigenvalue problems in coarse spaces 
Vhm (^ = 2, 3, ■ ■ ■ , n — 1). This is better than solving the eigenvalue problem on 
the finest finite element space directly, because solving source problem needs much 
less computation than solving the corresponding eigenvalue problem. 

Way 2. ("Multi-space method"): In this case, T4j, is defined on the same mesh 
Th but using higher order finite element than Vhf._^- In order to describe the scheme 
simply, we suppose the exact eigenfunction has sufficient regularity. We use the 
linear finite element space to solve the original eigenvalue problem (I2.12p on Vh, 
and solve the source problem (13. 4 p in higher order finite element space with the way 
that the order of Vh^ is one order higher than Vh^_^. Then we have the following 
error estimate for the final eigenpair approximation {\h^,Uh„) 

IA-A.J < r^a{Hf--'5l{\) = 0{H'-), (5.3) 

h-uuM < Va{Hr-'SH{X) = 0{H-). (5.4) 

The improved error estimates above just cost solving the source problems on the 
same mesh but in higher order finite element spaces and eigenvalue problems in the 
lowest order finite element space. 

6 Numerical results 

In this section, we give two numerical examples to illustrate the efficiency of the 
multi-level correction algorithm proposed in this paper. We solve the model eigen- 
value problem (II. ip on the unit square Q = (0, 1) x (0, 1). 

12 



6.1 Multi-space way 

Here we give the numerical results of the multi-level correction scheme in which the 
finer finite element spaces are constructed by improving the finite element orders on 
the same mesh. We first solve the eigenvalue problem fll.ll) in linear finite element 
space on the mesh Th- Then do the first correction step with quadratic element and 
cubic element for the second one. 

Here, we adopt the meshes which are produced by regular refinement from the 
initial mesh generated by Delaunay method to investigate the convergence behaviors. 
Figure [1] shows the initial mesh. Figure |2] gives the corresponding numerical results 
for the first eigenvalue Ai = 27r^ and the corresponding eigenf unction. From Figure 




Figure 1: Initial mesh for multi-space way 

[21 we can find each correction step can improve the convergence order by two for 
the eigenvalue approximations and one for the eigenfunction approximations with 
multi-space way when the exact eigenfunction is smooth. 

6.2 Multi-grid way 

Here we give the numerical results of the multi-level correction scheme where the 
finer finite element spaces are constructed by refining the existed mesh. We first 
solve the eigenvalue problem (11. ip by linear finite element space on the mesh Th- 
Then refine the mesh by the regular way such that the size of the resultant mesh 
hk = 0{H^) to obtain the mesh Th^, {k = 2,- ■ ■ ,n) and solve the auxiliary source 
problem f l3.4p in the finer linear finite element space T4j. defined on T/ij. and the 
corresponding eigenvalue problem in VH,hk- 

Figure E] gives the corresponding numerical results for the first eigenvalue A = 27r^ 
and the corresponding eigenfunction on the uniform meshes. From Figure [31 we 
can also find each correction step can improve the convergence order by two for the 
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Errors by Multiple Corrections: multi-space way 

-B-\-''\ 

— slope=-1 

— slope=-2 

— slope=-3 




Errors by Uultiple Corrections: multi-space way 

-B-Iivv^u|| 
-0-l|vu,'-vu|| 

— slope=-0.5 
slope=-1 
slope=-1.5 




10 10 

Number of elements 



10 10 

Number of elements 



Figure 2: The errors for the eigenvalue approximations by multi-level correction 
algorithm for the first eigenvalue 27r^ and the corresponding eigenfunction with 
multi-space way 

eigenvalue approximations and one for the eigenfunction approximations with the 
multi-grid way. 



7 Concluding remarks 

In this paper, we give a new type of multi-level correction scheme to improve the 
accuracy of the eigenpair approximations. We can use the better eigenvalue and 
eigenfunction approximation {Xh„,Uh„) to construct an a posteriori error estimator 
of the eigenpair approximation for the eigenvalue problem ([3 [T3]). 

Furthermore, our multi-level correction scheme can be coupled with the multigrid 
method to construct a type of multigrid and parallel method for eigenvalue problems 
(|22]). It can also be combined with the adaptive refinement technique for the 
singular eigenfunction cases. These will be our future work. 
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